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^ ■ Abstract 
(N 

' The usual theoretical treatments of the near-threshold pp — > pp-n^ reac- 

m 

\ tion are based on various phenomenological Lagrangians. In this work we 

<N . 

examine the relationship between these approaches and a systematic chiral 

in 

\ perturbation method. Our chiral perturbation calculation indicates that the 

pion rescattering term should be significantly enhanced as compared with 
^ ' the traditional phenomenological treatment, and that this term should have 

c ■ 

substantial energy and momentum dependence. An important consequence 



^ \ of this energy- momentum dependence is that, for a representative threshold 



kinematics and within the framework of our semiquantitative calculation, the 
rescattering term interferes destructively with the Born-term in sharp contrast 
to the constructive interference obtained in the conventional treatment. This 
destructive interference makes theoretical cross sections for pp ppn^ much 
smaller than the experimental values, a feature that suggests the importance 
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of the heavy-meson exchange contributions to explain the experimental data. 
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I. INTRODUCTION 



Recently Meyer et al. [0 carried out high-precision measurements of the total cross 



sections near threshold for the reaction 



p+p^p+p+n. 



(1) 



These measurements were confirmed by Bondar et al. |0. The early theoretical calculations 
underestimate these s-wave 7r° production cross sections by a factor of ~5. The basic 



features of these early calculations may be summarized as follows. The pion production 
reactions are assumed to be described by the single nucleon process (the Born term), Fig. 1(a), 
and the s-wave pion rescattering process. Fig. 1(b). The vr-A^ vertex for the Born term is 
assumed to be given by the pseudovector interaction Hamiltonian 



where gA is the axial coupling constant, and /tt = 93 MeV is the pion decay constant. 
The first term represents p-wave pion-nucleon coupling, while the second term accounts for 
the nucleon recoil effect and makes Tio "Galilean-invariant". For s-wave pion production 
only the second term contributes. Since this second term is smaller than the first term 
by a factor of ~ 171^^/771^-, the contribution of the Born term to s-wave pion production is 
intrinsically suppressed, and as a consequence the process becomes sensitive to two-body 
contributions. Fig. 1(b). The s-wave rescattering vertex in Fig. 1(b) is commonly calculated 
using the phenomenological Hamiltonian 



The two coupling constants Ai and A2 in Eq.(|^) were determined from the and 6*31 pion 
nucleon scattering lengths ai/2 and 03/2 as 





Til = 4:11 — —-ipTz-nip + An — "^ipT-nXTri/j 



(3) 




(4a) 
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A2 = ^ ( 1 + — ) (ai/2 - 03/2). (4b) 

The current algebra prediction for the scattering lengths, ai/2 = — 203/2 = mT^/ATrf^ = 
0.175m~^, implies that only chiral symmetry breaking terms will give a non- vanishing value 
of the coupling constant Ai in Eq.(^. Therefore Ai is expected to be very small. Indeed, the 
empirical values ai/2 — 0.175m^^ and 03/2 — — O.lOOm^^ obtained by Hohler et al. [0 lead to 
Ai ~ 0.005 and A2 ~ 0.05. So the contribution of the Ai term in Eq.(|^) is significantly sup- 
pressed. Meanwhile, although A2 is much larger than Ai, the isospin structure of the A2 term 
is such that it cannot contribute to the 7r° production from two protons at the rescattering 
vertex in Fig. 1(b). Thus, the use of the phenomenological Hamiltonians, Eqs.(0) and (^, 
to calculate the Born term and the rescattering terms illustrated in Figs. 1(a) and (b), gives 
significantly suppressed cross sections for the pp ppir^ reaction near threshold. Therefore, 
theoretically calculated cross sections can be highly sensitive to any deviations from this 
conventional treatment. These delicate features should be kept in mind in discussing the 
large discrepancy (a factor of ~5) between the observed cross sections and the predictions 
of the earlier calculations. 

A plausible mechanism to increase the theoretical cross section was suggested by Lee 
and Riska They proposed to supplement the contribution of the pion-exchange diagram. 
Fig. 1(b), with the contributions of the short-range axial-charge exchange operators which 
were directly related to heavy- meson exchanges in the nucleon-nucleon interactions 0. Ac- 
cording to Lee and Riska, the shorter-range meson exchanges (scalar and vector exchange 
contributions) can enhance the cross section by a factor 3-5. Subsequently, Horowitz et al. 
p!0| demonstrated, for the Bonn meson exchange potential, a prominent role of the a me- 
son in enhancing the cross section, thereby basically confirming the conclusions of Lee and 
Riska. The possible importance of heavy-meson exchanges may be inferred from the follow- 
ing simple argument. Consider Fig. 1(b) in the center-of-mass (CM) system with the initial 
and final interactions turned off and with the exchanged particle allowed to be any particle 
(not necessarily a pion). At threshold, go = ^n, 9 = 0, p'l = p'2 = 0, so that any exchanged 
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particle must have ko = = 70 MeV and |fc| = ^Jm-^rriN -\- {'rnT^/2y ~ 370 MeV/c, 

which imphes k"^ = — 171^^171 n. Thus the rescattering process probes two-nucleon forces at 
distances ~ 0.5 fm corresponding to a typical effective exchanged mass y^m^J^rTJv = 370 
MeV. Its sensitivity to the intermediate-range A^-A^ forces indicates the possible importance 
of the two-body heavy meson axial exchange currents considered by Lee and Riska. The par- 
ticular kinematical situation we considered here shall be referred to as the typical threshold 
kinematics. 



Meanwhile, Hernandez and Oset [11] considered the off-shell dependence of the tiN s- 
wave isoscalar amplitude featuring in the rescattering process, Fig. 1(b). They pointed out 
that the s-wave amplitude could be appreciably enhanced for off-shell kinematics pertinent 
to the rescattering process. We have seen above that, for the typical threshold kinemat- 
ics, the exchanged pion can indeed be far off-shell. The actual kinematics of course may 
deviate from the typical threshold kinematics rather significantly due to energy-momentum 
exchanges between the two nucleons in the initial and final states, but the importance of 
the off-shell kinematics for the exchanged pion is likely to persist. Hernandez and Oset 
examined two types of off-shell extrapolation: (i) the Hamilton model for vrA^ isoscalar am- 



plitude based on cr-exchange plus a short range piece ||T^, and (ii) an extrapolation based 
on the current algebra constraints. In either case the enhancement of the total cross sec- 
tion due to the rescattering process was estimated to be strong enough to reproduce the 
experimental data. A more detailed momentum-space calculation carried out by Hanhart 
et al. [jl3[ supports the significant enhancement due to an off-shell effect in the rescattering 
process, although the enhancement is not large enough to explain the experimental data. It 
should be emphasized that Hanhart et a/.'s calculation eliminates many of the kinematical 
approximations employed in the previous calculations. 

Given these developments based on the phenomenological Lagrangians, we consider it 
important to examine the significance of these phenomenological Lagrangians in chiral per- 
turbation theory (xPT) [|l4|,p!5| which in general serves as a guiding principle for low-energy 
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hadron dynamics. In the present work we shall describe an attempt at relating the tradi- 
tional phenomenological approaches to xPT. The fact that xPT accounts for and improves 
the results of the current algebra also makes it a natural framework for studying thresh- 
old pion production. Furthermore, in this low-energy regime, it is natural to employ the 
heavy- fermion formalism (HFF) |T^. The HFF has an additional advantage of allowing easy 
comparison with Eqs.(|D and @. 

It should be mentioned, however, that the application of xPT to nuclei involves some 



subtlety. As emphasized by Weinberg naive chiral counting fails for a nucleus, which 
is a loosely bound many-body system. This is because purely nucleonic intermediate states 
occurring in a nucleus can have very low excitation energies, which spoils the ordinary chiral 
counting. To avoid this difficulty, one must first classify diagrams appearing in perturbation 
series into irreducible and reducible diagrams, according to whether or not a diagram is free 
from purely nucleonic intermediate states. Thus, in an irreducible diagram, every interme- 
diate state contains at least one meson. The xPT can be safely applied to the irreducible 
diagrams. The contribution of all the irreducible diagrams (up to a specified chiral order) is 
then to be used as an effective operator acting on the nucleonic Hilbert space. This second 
step allows us to incorporate the contributions of the reducible diagrams. We may refer 
to this two-step procedure as the nuclear chiral perturbation theory (nuclear xPT). This 
method was first applied by Weinberg [|l^] to chiral-perturbation-theoretical derivation of 



the nucleon-nucleon interactions and subsequently used by van Kolck et al. [|T8[ . Park, Min 
and Rho (PMR) |jl9[ applied the nuclear xPT to meson exchange currents in nuclei. The 
success of the nuclear xPT in describing the exchange currents for the electromagnetic and 
weak interactions is well known [p!9|-pl|. The present paper is in the spirit of the work of 
PMR. 

This article is organized as follows: In the next section we define our pion field and the 
chiral counting procedure. Then in section III we present the two lowest order Lagrangians, 
discuss their connection to the early works on this reaction and determine within certain 
approximations the numerical values of the effective pion rescattering vertex strength, Kth- 
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In section IV we briefly discuss the connection between the transition matrix for this reaction 
and the xPT calculated amphtude. In section V we present necessary loop corrections to 
the Born term, and in section VI we calculate the cross section and discuss the various 
approximations and the uncertainties of the low energy constants in xPT. Finally in section 
VII, after discussing some higher chiral order diagrams, we present our main conclusions. 
A work very similar in spirit to ours has recently been completed by Cohen et al. ||2^ . 



II. CHIRAL PERTURBATION THEORY 

The effective chiral Lagrangian £ch involves an SU(2) matrix U{x) that is non-linearly 
related to the pion field and that has standard chiral transformation properties An 



example is [24 



U{x) = ^1 - [k{x)/UY + ^r-7T{x)/U. (5) 

In the meson sector, the sum of chiral-invariant monomials constructed from U{x) and its 
derivatives constitutes the chiral- symmetric part of Cch- Furthermore, one can construct 
systematically the symmetry-breaking part of £ch with the use of a mass matrix M. the 
chiral transformation of which is dictated by that of the quark mass term in the QCD 
Lagrangian. To each term appearing in £ch one can assign a chiral order index v defined by 

v = d-2, (6) 

where d is the summed power of the derivative and the pion mass involved in this term. 
A low energy phenomenon is characterized by a generic pion momentum Q, which is small 
compared to the chiral scale A ~ 1 GeV. It can be shown that the contribution of a term of 
chiral order v carry a factor {Q/AY, where Q represents either Q or the pion mass m^r. This 
suggests the possibility of describing low-energy phenomena in terms of Cch that contains 
only a manageably limited number of terms of low chiral order. This is the basic idea of 
XPT. 
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The heavy fermion formahsm (HFF) allows us to easily extend xPT to the meson- 
nucleon system. In HFF, the ordinary Dirac field ip describing the nucleon, is replaced by 
the heavy nucleon field N{x) and the accompanying "small component field" n{x) through 
the transformation 

ip{x) = exp{—imNV ■ x) [N{x) + n{x)] (7) 

with 

^iV = N, fn = -n, (8) 



where the four-velocity t>^ is assumed to be almost static, i.e., ~ (1,0,0,0) p5|. Elim- 
ination of n{x) in favor of N{x) leads to expansion in dfj_/mi\i. Since m^v ~ 1 GeV ~ A, 
an expansion in dfj_/mi\f may be treated like an expansion in d^/A. Cch in HFF consists 
of chiral symmetric monomials constructed from U{x), N{x) and their derivatives and of 
symmetry-breaking terms involving Ai. The chiral order u in HFF is defined by 

= d + n/2 -2, (9) 

where d is, as before, the summed power of the derivative and the pion mass, while n is the 
number of nucleon fields involved in a given term. As before, a term in Cch with chiral order 
u can be shown to carry a factor {Q/AY ^ 1. In what follows, u stands for the chiral order 
defined in Eq.@. 

In addition to the chiral order index u defined for each term in Cch, we assign a chiral 
order index u for each irreducible Feynman diagram appearing in the chiral perturbation 
series for a multifermion system [|1^]. Its definition is 



4-^jv-2C + 2L + ^z/i, (10) 



where Ejsf is the number of nucleons in the Feynman diagram, L the number of loops, and C 
the number of disconnected parts of the diagram. The sum over i runs over all the vertices 
in the Feynman graph, and Ui is the chiral order of each vertex. One can show that an 
irreducible diagram of chiral order u carries a factor {Q/AY <^ 1. 
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In the literature the term "effective Lagrangian" (or "effective Hamiltonian" ) is often used 
to imply that that Lagrangian (or Hamiltonian) is only meant for calculating tree diagrams. 
The Hamiltonians given in Eqs.(|^) and (^ are regarded as effective Hamiltonians of this 
type. We must note, however, that the effective Lagrangian in xPT has a different meaning. 
Not only can £ch be used beyond tree approximation but, in fact, a consistent chiral counting 
even demands inclusion of every loop diagram whose chiral order v is lower than or equal 
to the chiral order of interest. As will be discussed below, for a consistent xPT treatment 
of the problem at hand, we therefore need to consider loop corrections. However, since the 
inclusion of the loop corrections is rather technical, we find it useful to first concentrate on 
the tree-diagram contributions. This simplification allows us to understand the basic aspects 
of the relation between the contributions from xPT and the phenomenological Hamiltonians, 
Eqs.(0) and (|^). Therefore, in the next two sections (HI and IV) we limit our discussion to 
tree diagrams. A more elaborate treatment including loop corrections will be described in 
section V. 



III. TREE DIAGRAM CONSIDERATIONS 



In order to produce the one-body and two-body diagrams depicted in Figs. 1(a) and 1(b), 
we minimally need (see below) terms with u = 1 and 2 in Cch- We therefore work with 



(11) 



where C^'^'^ represents terms of chiral order u. Their explicit forms are [p!5|,p€ 



£(0) _ flTrld^Wd^U + ml{W + U - 2)] 
+ N{iv -D + g^S- u)N 



(12a) 
(12b) 
(12c) 



^9a 



2m^-N{S-D,v.u}N 
+ 2cimlNNTT{U + W - 2) 



(12d) 
(12e) 



+ {c2-^)N{vufN + c^Nu-uN (12f) 
-^iNN){mS-uN) (12g) 
- ^{NS^'N){Niu^N) (12h) 



2mjv 



In the above 

i ^ ^Juix), (13) 

u,^i{S}d,i-id,e), (14) 

D^N^id. + Wi'^^d.m. (15) 
and 5^ is the covariant spin operator defined by 

^,^i75[?^,7M]. (16) 

In C^^^ above we have retained only terms of direct relevance for our discussion. The coupling 
constants ci,C2 and C3 can be fixed from phenomenology [15]. They are related to the 
pion-nucleon cr-term, a^Nit) ~ {p'\Th{uu + dd)\p) {m = average mass of the light quarks, 
t = ip' — p)^), the axial polarizability a a and the isospin-even nN s-wave scattering length 
a+ = |(ai/2 + 203/2) ~ — O.OOSm^^ (The explicit expressions will be given below.) It 
should be noted that in HFF, a part of the term in C^^^ with the coefficient (c2 — g^/Sm^), 
namely the —g^/SniN piece, represents the s-wave n-N scattering contribution, which in a 
traditional calculation is obtained from the crossed Born-term . 

The four- Fermi non- derivative contact terms in Eq.(|l^) were introduced by Weinberg 



1% and further investigated in two- and three-nucleon systems by van Kolck et al. [[18 



Although these terms are important in the chiral perturbative derivation of the nucleon- 



nucleon interactions ||r^|18|, they do not play a major role in the following discussion of the 
threshold pp —>■ ppn^ reaction. We therefore temporarily ignore these four-fermion terms 
and come back to a discussion of these terms in the last section. 
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The Lagrangian (0) leads to the pion-nucleon interaction Hamiltonian 



where 



^lS = |fiV[o-V(T-7r)]iV + -i^iVT-7rx^iV (18a) 



+4[2cim27r2-(c2--^)7r2-C3(97r)2]iViV 



;i8b) 



Here Tij-^] represents the term of chiral order u. 



We now compare Ti-int resulting from xPT, Eq.([T7|) with the phenomenological effective 
Hamiltonian TYq + "Hi, Eqs.(|]) and (0). (The reader is reminded that the chiral index u 
should not be confused with the suffix appearing in Tio and Tii.) Regarding the uNN 
vertices, we note that the first term in ?-^*^°^ and the first term in Ti^^^ exactly correspond 
to the first and second terms, respectively, in TYq- Thus the so-called Galilean-invariance 
term naturally arises as a l/mj^ correction term in HFF. As for the tt-rNN vertices, we can 
associate the second term in 7ij°j to the A2 term in Hi, and second term in to the Ai 
term in Hi. This suggests the following identifications: 

= 7^ (19) 

and 

ArrXi/m^ = f [2c, - (c, - ^g^)^ - Cag] ^ .{k,g) (20) 

In Eg. (pOD , q = (cj^, q) and k = (ti;fc,k) stand for the four-momenta of the exchanged- 
and final pions, respectively, see Fig. 1(b). Since, as already discussed, the A2 term is not 
important for our purposes, we shall concentrate on the Ai term. The best available estimates 



of the coefficients q (i=l-3) can be found in Refs. |[T5|j27|| , which give 
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Cl 



4m?. 



-0.87 ± 0.11 GeV^\ 



(21a) 
(21b) 



C3 = TT 



2 



«A + ^(S + .i) 

J TT 

5.25 ± 0.22 GeV-\ 



(21c) 
(21d) 



C2 



J TT 

27712- 



47r(l + ^)a^ 

+ 2Ci - C3 + 



647r/4 



3.34 ± 0.27 GeV 



-1 



(21e) 
(21f) 
(21g) 



The numerical results are based on the experimental values: o"7rAr(0) = 45 ± 8 MeV ^8 



ua = 2.28 ± 0.10m^3 0, and a+ = (-0.83 ± 0.38) ■ lO-^m^^ [||]. We shall show in section 
VI that the uncertainties in the numerical value for C2 might be larger than quoted in 
Eq.( |21^ ). In fact, the terms in Eqs. (|21b| )- (|21g|) proportional to the {qa/ f-wY are (9((m,r/A)^) 
corrections arising from finite terms of C^'^\ However, since the present section is just an 
introduction to a later systematic treatment, this inconsistency in "accuracy" will be ignored 
for the moment. 

Now, for on-shelllow energy pion-nucleon scattering, i.e., fc~g~(m7r,0), we equate 



AnXi/m.,, = = K{k = {ttIt, , 0) , q = {mj, , 0)) , 



(22) 



where 



From Eq.(^) we have 



J TT 



J TT 



c + 



a: 



c = 2ci - C2 - C3, 



2ml 



47r(i + IlhLy 

rriN 



^9aK 
647r/4 



(23) 
(24) 



(25) 



which results in 
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«o = -27rri + ^V^ + M!|. (26) 

V m^y 1287r ^ > 

The above cited empirical value for a"*" leads to 

c = (0.59 ± 0.09) GeV"^ (27) 

Ko = (0.87 ± 0.20) GeV"\ (28) 

We now interpret these results in terms of Ai of Eq.(|^). Conventionally, Ai is determined 
from Eg. (|^) which is the first term in Eq.(^). Thus 



which gives 

47rAi 



'^^^ -2vr fl + ^ ) a\ (29) 



(0.43 ± 0.20) GeV~\ (30) 



or Ai = 0.005 ± 0.002. This is the "standard value" used in the literature [|7|,p0|. On the 
other hand, the r.h.s. of Eg . (|22D based on xPT gives from Eq. (^8D 

47rAi 



(0.87 ± 0.20) GeV^\ (31) 



which is about twice as large as the conventional value. This means the second term in 
Eq. ( p6D is almost as large as the first term. Thus xPT leads to a substantial modification of 
the commonly used formula, Eq. (^) or Eq. ( |29D . This large "higher chiral order" corrections 
due to C^"^^ [the term proportional to {qa/ f-wY' in Eq.(PBD] indicates that xPT does not 
converge very rapidly in this particular case. This apparent lack of convergence is probably 
due to the fact that the first terms in expansion, the vr-A^ isoscalar scattering length a"*", is 
exceptionally small. 

To develop further the connection between the traditional and the xPT approaches, we 
return to a discussion of Eq.(20). Obviously, the constant Ai cannot be fully identified 



with k(A;, q) which depends on the momenta q and k. In fact, the momentum dependence of 
k(/c, q) should play a significant role in describing the physical pion-nucleon elastic scattering 
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process where Uq = \Jml. + q^, uju = yrn^ + k^. An additional crucial point in the present 
context is that, in the rescattering diagram Fig. 1(b), the exchanged pion can be far off-shell, 
and therefore the q and k dependence in K,{k, q) may play an even more pronounced role. 
As an illustration, let us consider again the typical threshold kinematics discussed in the 
introduction: q ~ (m^, 0) and k ~ (Im^r, ^/rf^jn^). If we denote by Kth the value of K,{k, q) 
[Eq. (pOD ] corresponding to the typical threshold kinematics, we have 



2 

ml 



Hth — „2 

J TT 



1 / 9a\c3 



2ci - - C2 



2 \ SniNj 2 
The use of the central values for the coupling constants Ci, C2 and C3 leads to 



(32) 



AnXi/m^ = Kth ~ -1.5GeV~\ (33) 

Thus the strength of the s-wave pion-nucleon interaction here is much stronger than the 
on-shell cases, see Eqs.(pUD and (pi]), and the sign of the off-shell coupling strength is op- 
posite to the on-shell cases. The first feature is qualitatively in line with the observation of 
Hernandez and Oset ||ll| that the rescattering term should be larger than previously consid- 
ered. However, the sign of the typical off-shell coupling in our case [Eq.(|33|)] is opposite to 
the one used in Ref. [|ri|]. As will be discussed later, this flip of the sign drastically changes 
the pattern of interplay between the Born and rescattering terms. We must emphasize that 
the off-shell enhancement depends strongly on the values of Ci, C2 and C3, which, as discussed 



in Ref. [|T5|,^, are not known very accurately. It is therefore important to examine to what 
extent the existing large ambiguities in ci,C2 and C3 affect the off-shell enhancement of the 
pp ppTc^ reaction. We shall address this question in section VI. 



IV. TRANSITION OPERATORS FOR pp ppn^ 

As explained earlier, in the nuclear xPT we first use xPT to calculate the contributions 
of the irreducible diagrams. Let T represent the contributions of all irreducible diagrams 
(up to a specified chiral order u) for the pp —>■ pprr^ process. Then we use T as an effective 
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transition operator in the Hilbert space of nuclear wavefunctions. Consequently, the two- 
nucleon transition matrix element T for the pp — > ppn^ process is given by 

T = ($/|r|$,), (34) 

where (|$/)) is the initial (final) two-nucleon state distorted by the initial-state (final- 
state) interaction. These distorted waves should be obtained by solving the Schrodinger 
equation with nucleon-nucleon interactions generated by irreducible diagrams pertinent to 
nucleon-nucleon scattering, thereby incorporating an infinite number of "reducible" ladder 
diagrams. In this section we concentrate on the derivation of the transition operator T, 
relegating the discussion of T and Eq.(|3^) to section VI. 
We decompose T as 

T = ^r('^), (35) 

where T*^^) represents the contribution from Feynman diagrams of chiral order u, as defined 
in Eq.(|IUp. The lowest value of u occurs for the Born term shown in Fig.2(a). For s-wave 
vr production at threshold the NNn vertex with z/ = 0, the first term in Eq.(|T^), cannot 
contribute; hence the lowest u for NNn vertex involving an external pion must be z/ = 1. 
The first term in Eq. (|18b|) provides this vertex. According to Eq.(0), the chiral order of 
Fig. 2 (a) is given byz/ = 4 — 2 — 2-2 + 2- + 1 = —1. As can be checked easily, there 
are no diagrams with z/ = since in the rescattering diagram, Fig.2(b), the second term in 
Eq.( |18a| ), which gives the NNtttt vertex with z/ = 0, is not operative here due to the isospin 
selection rule. The rescattering diagram in Fig. 2(b) with the indicated value of z/ at each 
vertex contributes to T^'^^^'^ . It should be noted that because of the —2C term in the chiral 
counting expression, Eq. (p!0|) , exchange-current-type diagrams such as Fig. 2(b) give higher 
values of z/. In this work we truncate the calculation of the transition operator T at z/ = 1. 
Thus, 

T = T(-^) + T(^) (36) 
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The above enumeration is, as briefly discussed in section III, far from complete because 
loop diagrams and counter terms and finite terms from have been left out. In Fig.3 we 
show the loop corrections to the Born term [Fig.2(a)]. The diagrams in Fig.3 all have u = 1 
and hence are of the same chiral order as the leading order rescattering diagram. Fig. 2(b). 
As discussed earlier, for the pp — ^ ppir^ reaction at threshold the contribution of the Born 
term is numerically suppressed so that the rescattering diagram, which is formally of higher 
chiral order by two units of z/, plays an essential role. This implies that a meaningful and 
consistent xPT calculation of this reaction must include the loop corrections to the leading- 
order Born term. However, we continue to postpone the discussion of loop corrections to 
the next section. 

The tree diagrams contributing to Eg. (|36D , Figs. 2(a) and 2(b), are as follows. The Born 
term, Fig.2(a), contributes to T^^^^ and the rescattering term. Fig. 2(b), contributes to T^^\ 
These contributions are given, respectively, by 

^Res = _9A^ ^^^^^ Q)Tp^^f^ (37b) 

where Pi and p'i [i = 1, 2) denote the initial and final momenta of the z-th proton, ki = 
Pi — p'i'i and K{ki,q) is as defined in Eq. (^OD . 

V. LOOP DIAGRAMS 

We have emphasized above that the loop corrections to the Born diagram, Fig.2(a), which 
has chiral order u = —1, are of the same chiral order u = 1 as the two-body pion rescattering 
process. Fig. 2(b). These loop corrections therefore must be included in a consistent i/ = 1 
calculation. 

For our present purposes it is not necessary to go into a general discussion of the renor- 
malization of the parameters in Cch- Instead we concentrate on an estimation of the size of 
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the finite loop corrections to the specific tree level terms shown in Fig. 2. This will be done 
by applying standard Feynman rules and using dimensional regularization . Specifically, 
we only need consider the loop corrections to the single n^NN vertex in the s-wave channel: 



/T-Born 1 /T-Corr 

^-1 "T ^+1 



9a 



(38) 



where Si = (0, ^cTj) is the spin of the i-th proton and V is the amplitude to be calculated. 
For the Born term [Fig.2(a)] itself we have: 



Voa = 1 



(39) 



given by Eq. (|37a|) . The loop diagrams [Figs.3(a)-(f)], which renormalize the s-wave Born 
term, give the following contributions: 



3a 



4 \U. 



vq 



v 



3b 



V 



3c 



1 J2(V) - M'^P) 



2// 



vq 



3d 



2// 



3A^ + [{vp)Jo{vp) + {vp')Jo{vp')] 



(vpyJoivp) - {vp'YJoivp') 



Here we have adopted the notations of Ref. |15]. Thus 



A, 



i J (27r)^m^2 _/2 



vq 



(40a) 



(40b) 



(40c) 



(40d) 



2m^^ L 



In 



167r2 A 



(41) 
(42) 



where the divergence is included in 



-4 r 



167r2 



+ -{-fE - 1 - ln47r) 



(43) 



d-A 2 

In this expression A denotes the dimensional regularization scale and 'Je = 0.557215. Fur- 
thermore, Jo and J2 in Eqs. di^) are defined by 
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M.) = -41. + ^ (l - 21n - jLy^^TT^arocos (^) (44) 



and 



2 



^2)Jo(^)_c^aJ (45) 



d- 1 

The two contributions to V, Eqs. ( [40c]) and ( ^40d| ) , originate from two different combinations 
of terms in Eqs. ([T8|) . To calculate Eq.( [40c|) , the second term in Eq. (|18a|) and the first term 
in Eq.( |18b|) are used at the vertices, whereas Eq.( [40d|) is calculated using the first term of 
Eq.( |18a| ) and the second term of Eq. ( |18b|) . 



The standard renormalization consists in the following procedure: 

(1) The loop contributions to V are separated into a divergent part, which we take to be 
proportional to L of Eq.(^) and which contains a pole at c? = 4, and a finite part: 

V3 = VsT + VsK^"^*" (46) 

(2) Local counter terms, which are of the same chiral order as the loop diagrams, are added. 
In our case these counter terms must come from the Lagrangian C^^\ 

- Ji^.^^'^0,N (47) 

to give two-nucleon diagrams with u = 1. The unknown constants Di are then written as a 
sum of a finite and an infinite part 

Di = AK™*'(A) + {^TrfS^L. (48) 

The constants Si are determined by requiring that the infinite part of Di cancel the divergent 
part Vsl'^. The remaining finite contributions which should be added to the Born term via 
Eq.(p8D, are 

V/oop = V3K*"^*^ + V,,i.K^™*^ (49) 

The amplitude V3 contains energy-independent and energy-dependent parts, as can be seen 
in Eq. (|40|) . The energy-independent part can be absorbed in the renormalization of the 



following physical parameters: the pion wave function renormalization factor [Fig.3(e)], 
the nucleon mass mjv and the nucleon wave function renormalization factor [(Fig.3(f)], 
as well as the axial coupling constant qa [Figs. 3(a), (b),(e),(f)]. For the evaluation of the 
energy-dependent part we use the typical threshold kinematics: vq = m-Tr, vpi = m7r, vp2 = 
0. Putting these values into the corresponding terms in eq. (|ID|), we obtain as the total 
contribution of the diagrams in Fig. 3 

Thus Va]-^*"**^ amounts to 10% of the Born term [Fig.2(a)]. In addition we have finite 
contributions from the counterterms of C^^\ Vc.t.K*"**^. We note that only very few of the 
low energy constants in the counterterms Dj|'^*"**^(A) are known |T5|. Some of the low 
energy constants in C^'^\ i?i|-^*"**'^(A) have been estimated in Ref. |TH] assuming A resonance 



saturation. The result indicates i?j|-^*™*'^(A) ~ (9(0.1). For an estimate of the low energy 
constants AK'™*^(A) in £(3) 

seems reasonable to assume that they are of the same order of 
magnitude as the _Bi|-^™*'^(A) in To be conservative let us assume L'j|'^™*'^(A) ^ 0{1); 
then we expect Vct.^*"**'^ ~ 0.1. It is clear that, if those coefficients were "unreasonably 
large" , the convergence of the whole chiral series would be destroyed. 

Altogether, after renormalization the total contributions from the loop terms are ex- 
pected to amount to at most 20% of the Born term. This is not a completely negligible 
contribution in the present context because, as will be discussed in the next section, there 
can be a significant cancelation between the Born and the rescattering terms. Nevertheless, 
since our present treatment involves other larger uncertainties, we will neglect the renormal- 
ization of the Born term and henceforth concentrate on the bare Born term [Fig. 2 (a)] and 
the rescattering term [Fig.2(b)]. 

VI. CALCULATION OF THE TWO-NUCLEON TRANSITION MATRIX 



We derived in section IV the effective transition operator T arising from the tree diagrams 
and, in section V, we estimated the additional contributions due to the loop corrections and 
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presented an argument for ignoring the loop corrections in this work. These considerations 
lead to the the HFF expression of T up to order z/ = 1, given in Eqs.(^) and (0), and this 
T is to be used in Eq.(0) to obtain the two-nucleon transition matrix T. 

A formally "consistent" treatment of Eq.(p^) would consist in using for |<l>j) and 1$/) 
two-nucleon wave functions generated by irreducible diagrams of order up to z/ = 1. A 
problem in this "consistent" xPT approach is that the intermediate two-nucleon propagators 
in Fig.l can be significantly off-mass-shell, which creates a difficulty in any xPT calculation. 
Another more practical problem is that, if we include the initial- and final- two-nucleon {N- 
N) interactions in diagrams up to chiral order = 1, these A^-A^ interactions are not realistic 
enough to reproduce the known A^-A^ observables. A pragmatic remedy for these problems 
is to use a phenomenological A^-A^ potential to generate the distorted A^-A^ wavef unctions. 
Park, Min and Rho |21|] used this hybrid approach to study the exchange-current in the 
n + p ^ '-f + d reaction and at least, for the low- momentum transfer process studied in Ref. 
|2T| , the hybrid method is known to work extremely well. 

Apart from the above-mentioned problem, there is a delicate aspect in the derivation 
of an effective two-body operator from a given Feynman diagram. Ordinarily, one works 
with r-space transition operators acting on r-representation wavef unctions, for the nuclear 
wavefunctions are commonly given in this representation. To this end, a Feynman ampli- 
tude which is most conveniently given in momentum space, is Fourier-transformed into the 
r-representation. This method works best for low momentum transfer processes which have 
substantial transition amplitudes for on-shell initial and final plane-wave states. However, 
the pp — s> ppn'^ reaction at threshold does not belong to this category. For this reaction it 
is essential to recognize that the nucleon lines that appear as external lines in Fig. 2 are in 
fact internal lines in larger diagrams illustrated in Fig.l. These internal lines can be far 
off-shell due to the initial- and final-state interactions. Indeed without this off-shell kine- 
matics, the Born term [Fig.2(a)] would not contribute at all ! In the conventional approach, 
however, one ignores this feature in deriving T in coordinate representation. For example, 
in Fourier-transforming an operator of the type of T^{^^, Eq. (|37b|) , even though pi and p'i 
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in Eg. (|37b| ) in fact can be anything due to momentum transfers caused by the initial and 
final A^-A^ interactions, it is a common practice to keep the energy of the propagating pion 
fixed at the value determined by the asymptotic energies of the nucleons. Hanhart et al. 
||r3| made a critical study of the consequences of avoiding these kinematical approximations. 
They worked directly with the two-nucleon wavefunctions in momentum representation. In 
the present work we do not attempt at detailed momentum-space calculations and sim- 
ply use the "conventional" Fourier transform method. Because of this and a few other 
approximations adopted, the numerical work presented here is admittedly of exploratory 
nature. Nonetheless, as we shall show, our semiquantitative study of T based on the chiral- 
theoretically motived transition operator T provide some valuable insight into the dynamics 
of the threshold pp — > ppn^ reaction. 

Let us denote the contribution of Fig. 2(b) for plane- wave initial and final states by 
< p[, P2, q\T^^^\pi, P2 >■ We first calculate this matrix element for the typical threshold 
kinematics described earlier; for the meson variables, q = (m^,0) and k = {171^^/2, k) with 
|k| = y/ni^jfij^. Correspondingly, the coupling strength n{k,q) [Eq.( ^7bD ] is taken to be 
= — l.SGeV^^ [Eq.(P3D]. Subsequently, by liberating the momentum variables pi, p'^, 
P2, and P2 from the on-mass-shell conditions (pf = m%, . . .), we treat <p'i, P2, q\T^^^\pi, p2 > 
as a function of pi, p'^, p2, and pg. Let T{p'^, p'2] pi, P2) stand for this function. We still 
require momentum conservation at each vertex, which imposes the conditions pi + P2 = 
Pi + P2 + Q = 0) k = p'^ — Pi = P2 — P2- ^(Pi) P'2'1 Pi' P2) can be easily Fourier 
transformed to give T^^'^ in r representation. The simplified treatment described here, which 
is commonly used in the literature, shall be referred to as the fixed kinematics approximation. 

Now, in the fixed kinematics approximation, T [Eqs.(^6]), (p7|)] is translated into differ- 
ential operators acting on relative coordinate of the two-nucleon wavefunctions: 

r^r = §^ — S-V„ (51a) 
f^r = -'^^^tH^-rf'ir), (51b) 
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where the derivative operator with subscript r is to act on the relative coordinate r between 
two protons, and S = |(cri — 0-2) . The trivial isospin operator has been dropped. 
The Yukawa function /(r) = exp(— yu'r)/47rr is defined with the effective mass fi' = ^m^r. 
We reemphasize that the simple Yukawa form /(r) arises only when the fixed kinematics 
approximation just discussed is used. 

From this point on, our calculation of T follows exactly the traditional pattern described 
in the literature. Thus T is evaluated by inserting the transition operators, T]^™" and T^f^, 
Eq.(^T]), between the initial and final nuclear states 

0i(r) = (v^/pr)mi,o(r)e'^i-«(47r)V2|3,3p^^^ 

(52) 

0y(r) = (l/pV)no,o(r)e^^o.«(47r)i/2|3,i5^), 

where p and p' are the asymptotic relative three-momenta of the initial and final two-proton 
systems. The wavefunctions are normalized as u^^ ^-^^ sin(pr — ^ttL + 5l,j) with 6lj being 
the N-N scattering phase shifts. For simplicity, the Coulomb interactions between the two 
protons is ignored. (The Coulomb force is known to reduce the cross section up to 30%. 0].) 
The explicit expression for the transition amplitude at threshold is obtained as 

TiEf) = 47r(^7^/Am^m3/2)(jBorn ^ ^Res)^ (53) 

Here, Ef = Ep/ + /2m-,j is the kinetic energy of the final state, and 

J^r = hm ^ rdrr-'^ f # + ^1 ^, (54a) 



pp' Jo r \dr r I r 



= hm 2k,,^^ / druo,of'{r)u^,o. (54b) 

p'^O pp' Jo 



The total cross section is obtained by multiplying the absolute square of the transition ampli- 
tude (averaged over the initial spins and summed over the final spins) with the appropriate 
phase space factor p{Ef) and the fiux factor 1/v: 

27r 



jdp{Ej)\T{Ef)\\ (55) 
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For a rough estimation one may approximate the energy dependence of the transition matrix 



as [pll] 



where a is the scattering length of the A^A^ potentiaL Then the cross section can be simply 
expressed as 

Otot= J\ A J?l{E,) (57) 



where 



J|2 = |jBorn + jRes|2^ (58a) 



rEj JEf — Ep/ 1 / Ep> 

I{Ef) = / dE,,\l (58b) 

Under the approximation (p6|) , the energy dependence of the cross section is solely given 



by I{Ef), which incorporates the phase space and the final state interaction effect (in the 
Watson approximation [0). 

We have calculated the integrals J^™" and J^f^ for representative nuclear potentials: 
the Hamada- Johnston (HJ) potential and the Reid soft-core potential (RSC) |Q. The 
results are given in Table I, and the corresponding cross sections are presented in Table II. 
These results indicate that, for the nuclear potentials considered here, the value of |J| is 
much too small to reproduce the experimental cross section. If we define the discrepancy 
ratio R by 

R ^ aLT/a^t , (59) 

with 0"^/ taken from Ref. then R = 80 [R = 210) for the Hamada- Johnston (Reid 
soft-core) potential, and R happens to be almost constant for the whole range of Ef < 23 
MeV for which cr^o/ is known. Thus, although the off-shell behavior of the s-wave pion 
scattering amplitude derived from the chiral Lagrangian does enhance the contribution of 
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the rescattering process over the value reported in the hterature, the sign change that occurs 
in /€ as one goes from kq [Eq.(p2D] to Hth [Eq.(p^)] results in a significant cancelation between 
the Born term J^°™ ^ind the rescattering term J^f^, leading to the very small cross sections 
in Table II |^^. The drastic cancelation between J^°™ ^^^1 J^f^ found here means that 



the calculated cross sections are highly sensitive to the various approximations used in our 
calculation and also to the precise values of the constants ci,C2 and C3 of Eq. f^Tl) . We will 
discuss these two questions in the next two paragraphs. 

We adopted the threshold kinematics approximation and neglected the energy- 
momentum dependence in Eq . (^Of) and treated the vertices in Figs.l and 2 as fixed numbers, 
i.e., K{k,q) = = constant. In addition, although the loop corrections of chiral order u 
= 1, shown in Fig. 3 (a), automatically introduces energy- momentum dependent vertices, we 
ignored this feature. The fact that the kinematics of the reaction Eq.(|I|) requires highly 
off-shell vertices leads to the expectation that the vertex form factors can be very important 
and invalidate the threshold kinematics approximation leading to Eq. (|5T|) . In this connec- 
tion we note that a momentum-space calculation [|T^, which is free from this approximation, 
indicates that even a negative value of Ai could lead to the moderate enhancement of the 
cross section. 

The strong cancelation between the Born and rescattering terms also means that, even 
within the framework of the fixed kinematics approximation, the large errors that exist in 
the empirical value of a'^ and the Ci,C2 and C3 constants can infiuence the cross sections 
significantly. To assess this infiuence, we rewrite Eq.(p2D as 

The use of the experimental values for a+ and Ci quoted earlier leads to 

Kth = (-1.5±0.4)GeV~\ (61) 

With this uncertainty taken into account, the ratio R ranges from R = 25 to R = 2100 for 
the Hamada- Johnston potential, and from R = 50 to R = 3.4 x 10^ for the Reid soft-core 
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potential. To further examine the uncertainties in the C^^^ constants we remark that the 
value of C2 + C3 can be extracted from the known pion-nucleon effective range parameter ft"*". 
The low energy pion-nucleon scattering amplitude is expanded as: 

/+ = «+ + b+q^ + ■■■ (62) 

where q is the pion momentum and = (— 0.044 ±0.007)m^^ [|^. If we use C^^^ to calculate 
the s-wave pion-nucleon amplitude we find: 

i+= 1 fl + :!!^V7!^Vfe + C3-A)A. (63) 



7^ i + — \-r 1C2 + C3-- — ) — 

2n V rriN/ \ Jtt J ouin 
and then Eq. (|32|) leads to 



^ 2m 



- 7rm^2 f 1 + ^) b+ 
= (-2.7±0.6)GeV"^ (64) 

Since C3 is given directly by the experimental quantity a a [Eq. (|21d|) ], we consider Eq.(|63D 



as an alternative input to determine C2 in terms of fe"*" and C3. Then Eqs. (|63|) ,( 21d ) and the 
experimental value of b^ [0 give 

C2 = (4.5 ± 0.7)GeV-^ . (65) 



We note that this value is larger than the one given in Eq. ( |21g| ), indicating that the de- 
termination of C2 requires further studies. With the new value of K^th given in Eq.(|5^ 
we find that the discrepancy ratio R [Eq.(^^] can be as small as ~ 10. ( In this case 
|jRes| ^ |jBorn|. ^y^q exact Cancelation between the Born and the pion rescattering term 
occurs for Kth ~ — 2GeV^^.) 

Without attaching any significance to the detailed numbers above, we still learn the 
extreme sensitivity of af"/^ to the input parameters and that, despite this high sensitivity, 
CTj^t'^ still falls far short of alof (within the framework of the fixed kinematics approximation) . 
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VII. DISCUSSION AND CONCLUSIONS 



In this work we have used xPT to calculate the effective pion-exchange current contribu- 
tion to the pp — > ppiT^ reaction at threshold. As stated repeatedly, our aim here is to carry 
out a systematic treatment of T up to chiral order u = 1 [see Eq.(|36|)]. However, in order 
to make contact with the expressions appearing in the literature let us consider a very 
limited number of u = 2 diagrams. To be specific, we consider a diagram in Fig. 2(b) but 
with the u = (p-wave) ttNN vertex replaced with a, v = 1 (s-wave) vertex. Then, instead 
of Eq.( |51b| ), we will obtain 



^Res _ 9a ^ 
■'■1+2 — ~~r'^th^' 

Jtt 



f/'(r) (2 + ^)+/(r)^V, 



(66) 



which is the two-body transition operator used in Ref. |Q. Thus, we do recover the usual 
phenomenological parameterization in xPT, but this is just one of many v = 2 diagrams. 
Our systematic u = 1 calculation excludes all u = 2 diagrams. 

We have also ignored the exchange current contributions from scalar and vector two- 



nucleon exchanges. Following the xPT of Refs. [0,0 the vector meson exchange is largely 
accounted for via the four-nucleon contact terms illustrated in Fig. 4(a). If we had retained 
the last two terms of Eg. (p^ , the pion-nucleon interaction Tij^^, Eq.( |18b| ), would have had 
an additional piece 

= -^iNN)iNa.Vir-7z)N) 

+ ^^V (iVerA^) ■ {NV{t-7t)N). (67) 

The Tiint term of Fig. 4(a) has a <t ■ q structure, which means it describes p-wave pion 
production and therefore does not contribute to the threshold pp —>■ ppir^ reaction. The 
s-wave pion production contact term, also belonging to the type of diagram illustrated 
in Fig.4(a), enters as a — recoil correction to and therefore is of chiral order u = 

2. Formally, the chiral order u = 2 diagrams have no place in the present calculation 
limited to u = 1. However, in view of the great current interest in the possible large 
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contribution of the heavy-meson exchange diagrams, we make a few remarks on the s-wave 
u = 2 contact terms depicted in Fig.4(a). We note that the coordinate representation of 
this contact term contains 6^{r). Meanwhile, in the threshold pp — > ppn^ reaction the initial 
two-nucleon relative motion must be in p-wave (because of parity) and so its wavefunction 
vanishes at r = 0. Thus, even in a chiral order u = 2 calculation, the contact term Fig.4(a) 
corresponding to s-wave pion production will play no role. Including meson loops corrections 
to these contact terms [an example illustrated in Fig. 4(b)] would smear out the (5-function 
behavior, allowing them to have a finite contribution to the threshold pp — > ppn^ reaction. 
This involves, however, diagrams of even higher chiral order than u = 2. Thus, in order 
to include the strong effective isoscalar-vector repulsion of the A^-A^ forces {u exchange) 



contained in the four-nucleon contact terms of Weinberg's and van Kolck's et a/.'s |jn 
xPT description, we have to go to chiral order u = 3. 

Meanwhile, one may picture the "effective heavy mesons" as generated by multi-pion 
exchange diagrams like those illustrated in Fig. 5. These diagrams, which necessarily contain 
loops, represent a very limited class of z/ > 3 diagrams. For example, an important part of the 
effective scalar exchange between two nucleons involve intermediate vr-vr s-wave interaction 
which requires at least two loop diagrams like Fig.5(c). Thus, if we are to interpret the 
heavy-meson exchange diagrams of Lee and Riska |^ in the framework of nuclear xPT, 
we must deal with terms with chiral order z/ > 3, which at present is beyond practical 
calculations. 

We now recapitulate the main points of this article. 

1. Using xPT in a systematic fashion we have shown that the contribution of the pion 
rescattering term can be much larger than obtained in the traditional phenomenolog- 
ical calculations. This fact itself supports the suggestion of Hernandez and Oset |TT 



that the off-shell s-wave pion-nucleon scattering should enhance the rescattering con- 
tribution significantly. However, the sign of the enhanced rescattering vertex obtained 
in xPT is opposite to that used in Ref . [O] , at least for the typical threshold kinematics 
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defined in the text. This sign change in the couphng constant nth leads to a destructive 
interference between the Born and rescattering terms instead of the constructive in- 
terference found in Ref. The significant cancelation between these terms give rise 
to the very small cross section for the near-threshold pp —>■ ppn^ reaction calculated 
in this work. Although our particular numerical results were obtained in what we 
call the fixed kinematics approximation, these results at least indicate that the large 



enhancement of a^^^ obtained in Ref. [|Il|] is open to more detailed examinations. 



2. The fixed kinematics approximation (which is commonly used in the literature) should 
be avoided. There are at least two reasons why this is not a good approximation 
for this reaction: (i) The initial- and final-state interactions play an essential role 
in the near-threshold pp — ppvr" reaction; (ii) The theoretical cross section within 
the framework of the Born plus rescattering terms is likely to depend on the delicate 
cancelation between these two terms. In a momentum space calculation [|I^, we can 
easily avoid the fixed kinematics approximation. Such a calculation will allow us to 
work with full off-shell kinematics, to incorporate the xPT form factors in the Born 
term, and to reduce ambiguities in our calculation down to the level of uncertainties 
in the input parameters in xPT and the chiral counter terms. 

3. Several works P, |T0i|l3| indicate that the two-nucleon scalar (sigma) exchange can be 



very important. We gave in the introduction a simple kinematical argument for its 
plausibility, and our dynamical calculation (albeit of semiquantitative nature) seems 
to indicate the necessity of the sigma exchange contribution in order to explain the 
observed cross sections for the threshold pp — * pprc^ reaction. It is of great interest to 
see to what extent an improved xPT calculation based on momentum-space represen- 
tation helps sharpen the conclusion on the necessity of the sigma exchange diagram. 
Such a calculation is now in progress. If it is established that the heavy meson ex- 
change diagrams play an essential role in the threshold pp ppir^ reaction, it seems 
that we must resort to a modified version of xPT, for a brute force extension of our 
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treatment to z/ > 2 seems extremely difficult. An attempt to include vector meson de- 
grees of freedoms explicitly can be found e.g. in Ref. |19[. A purely phenomenological 
approach as used in may also be a useful alternative. 
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FIGURES 

FIG. 1. Single nucleon process (Born term) [Fig. 1(a)] and pion rescattering process for the 
pp ppiT^ reaction near threshold. 2T+i,25+i^^ denotes the isospin and angular momenta of the 
initial and final states. 

FIG. 2. Tree graphs: The Born term [Fig. 2(a)] {v = —1) and the pion rescattering term 
[Fig.2(b)] (i. = l). 

FIG. 3. Loop corrections to the Born term. 

FIG. 4. Generic four-fermion-pion vertex (contact term) [Fig. 4(a)] and an example of a loop 
correction to a contact term [Fig.4(b)]. 

FIG. 5. A few higher order diagrams contributing to the effective two-nucleon scalar exchange 
in nuclear xPT. 
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TABLES 

TABLE L J Born ^^^^ jRes f^^. ^j^g threshold kinematics [Eqs.(^),(|54b|)], calculated with 
the Hamada-Johnsotn (HJ) and Reid soft-core (RSC) potentials. 





HJ 


RSC 


rBorn 


-0.672 


-0.515 


jRes 
•^+1 


+0.505 


+0.413 
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TABLE II. The total cross sections (in 


/x6) as functions of 77 = 




calculated with the 


Hamada-Johnston (HJ) and Reid soft-core 


(RSC) potentials. 






V 






(^RSC 


.03 


.0000 




.0000 


.06 


.0003 




.0001 


.09 


.0011 




.0004 


.12 


.0024 




.0009 


.15 


.0043 




.0016 


.18 


.0069 




.0026 


.21 


.0100 




.0037 


.24 


.0138 




.0052 


.27 


.0182 




.0068 


.30 


.0232 




.0087 


.33 


.0289 




.0108 


.36 


.0352 




.0131 


.39 


.0421 




.0157 


.42 


.0496 




.0185 


.45 


.0577 




.0215 


.48 


.0665 




.0248 


.51 


.0759 




.0283 


.54 


.0859 




.0320 


.57 


.0965 




.0360 


.60 


.1078 




.0402 
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